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General expressions for the longitudinal and transverse resistivities of single-crystalline cubic and 
tetragonal ferromagnets are derived from a series expansion of the resistivity tensor with respect 
to the magnetization orientation. They are applied to strained (Ga,Mn)As films, grown on (001)- 
and (113)A-oriented GaAs substrates, where the resistivities are theoretically and experimentally 
studied for magnetic fields rotated within various planes parallel and perpendicular to the sample 
surface. We are able to model the measured angular dependences of the resistivities within the 
framework of a single ferromagnetic domain, calculating the field-dependent orientation of the mag- 
netization by numerically minimizing the free-enthalpy density. Angle-dependent magnetotransport 
measurements are shown to be a powerful tool for probing both anisotropic magnetoresistance and 
magnetic anisotropy. The anisotropy parameters of the (Ga,Mn)As films inferred from the mag- 
netotransport measurements agree with those obtained by ferromagnetic resonance measurements 
within a factor of two. 
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I. INTRODUCTION 



Realization of ferromagnetism in HI-V semiconduc- 
tors by introducing high concentrations of magnetic ele- 
ments has motivated intense research on the dilute mag- 
netic semiconductor (Ga,Mn)As. This system is con- 
sidered a potential candidate for spintronic applications 
due to its compatibility with conventional semiconduc- 
tor technologyi^ In (Ga,Mn)As, magnetic Mn acceptors 
are predominantly incorporated on cation sites as Mn 2+ 
ions having a total spin of S — 5/2. (Ga,Mn)As is para- 
magnetic at room temperature and undergoes a transi- 
tion to the ferromagnetic phase at the Curie temperature 
Tc, where maximum values of up to ^170 K have been 
reported so far^ The ferromagnetism has been success- 
fully explained within the Zener mean-field model by an 
indirect Mn-Mn exchange interaction mediated by delo- 
calized holes* 3 - (Ga,Mn)As is grown by low-temperature 
molecular-beam epitaxy and, if necessary, subjected to 
post-growth annealing to reduce the density of compen- 
sating defects. Considerable progress has been made in 
understanding its structural, electronic, and magnetic 
properties. In particular, anisotropic magnetoresistance 
(AMR), 4 i 5 i 6 i 7 i 8 i 9 planar Hall effect (PHE)^. and mag- 
netic anisotropy (MA), 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 have been identi- 
fied as characteristic features, making (Ga,Mn)As poten- 
tially suitable for field-sensitive devices and non-volatile 
memories. These properties have been shown to be 



governed by several parameters, such as Mn concen- 
tration, hole density, strain, or temperature. Most of 
the work carried out on AMR in (Ga,Mn)As, however, 
has been restricted to special cases where the magnetic 
field was applied parallel or perpendicular to the layer 
and equations for the angular dependence of the longi- 
tudinal and transverse resistivities, describing the AMR 
and the PHE, respectively, have been given only for in- 
plane configuration and polycrystalline films. A com- 
prehensive theoretical model describing the resistivities 
as a function of arbitrary field orientation is still miss- 
ing. Moreover, (Ga,Mn)As layers are usually grown on 
GaAs(OOl) substrates and only little is known about 
the magnetic properties of films grown on high-index 
substrates, 16 ! 17 ! 18 ! 19 . 20 ! 21 

In this work, the longitudinal and transverse resistiv- 
ities of (Ga,Mn)As layers, grown on (001)- and (113)A- 
oriented GaAs substrates, are studied for arbitrarily ori- 
entated magnetic fields. The anisotropy of the resistivi- 
ties and the MA are experimentally probed by rotating 
the magnetic field H at fixed field strengths within differ- 
ent planes parallel and perpendicular to the sample sur- 
face. General expressions for the resistivities, holding for 
single-crystalline cubic and tetragonal ferromagnets, are 
derived from a series expansion of the resistivity tensor 
with respect to the direction cosines of the magnetization 
M. The measured data are well modeled by applying 
the expressions to the given experimental configurations, 
assuming the (Ga,Mn)As films to consist of a single fer- 
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FIG. 1: Orientation of the (001) and (113)A samples with 
respect to the crystallographic axes. 

romagnetic domain. Analytical expressions, widely used 
in the literature to describe the angular dependence of 
AMR and PHE, are shown to be inapproriate to single- 
crystalline materials. Finally, anisotropy parameters are 
estimated from the low-field magnetotransport data and 
compared with those obtained from ferromagnetic reso- 
nance (FMR) spectroscopy. 



II. EXPERIMENTAL DETAILS 

40-nm-thick (001) and (113)A (Ga,Mn)As films with 
Mn concentrations of ~5% were simultaneously grown 
by low-temperature molecular-beam epitaxy (MBE) in a 
RIBER 32 MBE machine on semi-insulating GaAs(OOl) 
and GaAs(113)A substrates mounted together on the 
same Mo holder. A conventional Knudsen cell and a 
hot-lip effusion cell were used to provide the Ga and Mn 
fluxes, respectively. A valved arsenic cracker cell was op- 
erated in the non-cracking mode to supply AS4 with a 
maximum V/III flux ratio of about 5. First, a 30- nm- 
thick GaAs buffer layer was grown at a temperature of 
T s ss 580° C (conventional growth temperature for GaAs), 
then the growth was interrupted and T s was lowered to 
~250°C. The Mn concentrations in the (Ga,Mn)As films 
were determined by flux measurements. 

For the magnetotransport and FMR studies, the (001) 
and (113) A samples were cleaved into small rectangular 
pieces with edges along [110] and [110], and along [332] 
and [110], respectively. The [110] and [110] directions of 
the (001) sample were identified by selective wet chemical 
etching and the orientation of the [001] crystal axis in the 
(113)A sample was determined by x-ray diffraction. The 
crystallographic orientations of the samples are shown in 
Fig. 2] Hall bars with Ti-AuPt-Au contacts were pre- 
pared on several pieces of the cleaved (001) and (113)A 
samples with the current direction along [110] and [332], 
respectively. The width of the Hall bars was 0.3 mm and 
the longitudinal voltage probes were separated by 1 mm. 
Hole densities of 3xl0 20 cm~ 3 for the (001) sample and 
1.2xl0 20 cm~ 3 for the (113)A sample were determined 
by means of high-field magnetotransport measurements 
(up to 14.5 T) at 4.2 K using an Oxford SMD 10/15/9 VS 
liquid helium cryostat with superconducting coils. Least 



squares fits were performed to separate the contributions 
of the normal and anomalous Hall effect. Curie tem- 
peratures of Tc ~ 65 K and 54 K, respectively, were 
estimated from the peak positions of the temperature- 
dependent sheet resistivities at 10 mT. For the angle- 
dependent magnetotransport measurements carried out 
at 4.2 K, the Hall bars were mounted on the sample 
holder of a liquid-He-bath cryostat, which was positioned 
between the poles of a LakeShore electromagnet system 
providing a maximum field strength of 0.7 T. The sam- 
ple holder possesses two perpendicular axes of rotation, 
allowing for an arbitrary alignment of the Hall bars with 
respect to the applied magnetic field H. Using a dc cur- 
rent density of 8x 10 2 Acm -2 , the longitudinal and trans- 
verse resistivities piong and ptrans were measured at fixed 
magnitudes /lo-ff = 0.1, 0.25 and 0.7 T of H while ro- 
tating its orientation. Prior to each angular scan, the 
magnetization M was put into a clearly defined initial 
state by raising the field to 0.7 T where M is nearly sat- 
urated and aligned with the external field. The field was 
then lowered to one of the above mentioned magnitudes 
and the scan was started. 

The FMR measurements were carried out at 5 K in a 
commercial Bruker ESP 300 electron paramagnetic res- 
onance spectrometer operated at a fixed frequency of 
ujhf/Ztt ~ 9.3 GHz (X-band). The spectrometer consists 
of a microwave bridge for the high-frequency radiation 
and an electromagnet providing a variable dc magnetic 
induction up to 1 T. To increase the sensitivity, lock-in 
techniques were used in which the dc induction is super- 
imposed by a 100 kHz modulation field of 3.2 mT. 



III. THEORETICAL OVERVIEW 

In our theoretical considerations the total mag- 
netic moment arising from the Mn-ion/holc spin com- 
plex is treated within the framework of the Stoner- 
Wohlfarth modelj 2 ^ i.e., for temperatures below Tc the 
whole (Ga,Mn)As layer is assumed to consist of a sin- 
gle homogeneous ferromagnetic domain. This simple 
model has been astoundingly successful in describing 
a large variety of magnetization-related phenomena in 
(Ga,Mn)As. Under the given experimental conditions de- 
scribed above, domain nucleation and expansion, which 
have been shown to accompany in-plane and perpendic- 
ular magnetization-reversal processes^i*^ are expected 
to play only a minor role. Accordingly, we may write 
the magnetization as a vector M = Mm where M de- 
notes its magnitude and the unit vector m its direction. 
In terms of the polar and azimuth angles 8 and tp, re- 
spectively, which are defined in Fig. ^ the components 
of m read as m x = sin 8 cos tp, m y = sin 8 sin tp, and 
m z = cos 8. The equations used in the discussion of the 
angle-dependent magnetotransport data can be written 
in a concise way by introducing the unit vectors j, n, 
and t, which specify the current direction, the surface 
normal, and an in-plane vector defined by t = n X j, re- 
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spectively. Throughout this work, all vector components 
refer to the cubic coordinate system with the [100], [010], 
and [001] directions of the crystal denoted by x, y, and 
z, respectively. 



Insertion of Eq. into Eqs. QJ and elementary vector 
algebra yields the general expressions 



Plo 



A + C(j ■ mf + {B-C)Y, 



■2 2 



A. Longitudinal and transverse resistivities 

In standard magnetotransport measurements the lon- 
gitudinal and transverse voltages, measured along and 
across the current direction, arise from the components 
E[ ong — j ■ E and E trans — t ■ E of the electric field E, 
respectively. Starting from Ohm's law E = p - J , where p 
represents the resistivity tensor and J = Jj the current 
density, the corresponding longitudinal resistivity pi on g 
(sheet resistivity) and transverse resistivity ptrans (Hall 
resistivity) can be written as 
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(1) 



In (Ga,Mn)As, as in many other ferromagnets, the resis- 
tivity tensor sensitively depends on the orientation of M 
with respect to the crystallographic axes^ Thus, in or- 
der to quantitatively model the measured resistivities in 
the general case of an arbitrarily oriented magnetization, 
a universal mathematical relationship between piong and 
Ptrans and the direction cosines m, of M has to be de- 
rived. For these purposes, we follow the ansatz of Birss 23 
and Muduli et ali^i and write the resistivity tensor p as 
a series expansion in powers of nij using the Einstein 
summation convention: 



Pij = a ij + a ki]m k + akiij-nikmi + 



(2) 



For cubic symmetry Td, most of the components , akij , 
of the galvanomagnetic tensors vanish and, consider- 
ing terms up to the second order, we obtain 



1 0\ / ml 
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(3) 
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with the resistivity parameters 

A = an + flii22 , B — aim — 01122 , 



C = a 232 3 , D = a\ 



23- 



(4) 



P™an S c = C(j ■m)(t-m) + (B-C)Y^ U h m 2 t 

i 

- D(n-m), (5) 

which apply to single crystalline ferromagnetic materials 
of cubic symmetry. The transverse resistivity in Eqs. (JSJ 
includes the contribution of the anomalous Hall effect, 
which correlates with the perpendicular component of 
M , but it does not account for the ordinary Hall effect. 
For magnetic field strengths pqH < 1 T and hole concen- 
trations p > 10 20 cm -3 as in our experiments, however, 
the maximum contribution of the ordinary Hall effect 
is poH/ep f» 6- 10~ 6 57 cm (e denotes the elementary 
charge) , and thus about two orders of magnitude smaller 
than the measured peak values of Ptrans ( see Section IV) . 

In the following, Eqs. © are applied to three different 
experimental configurations using the relation 



mi = ■ m) + U(t ■ m) + n^n ■ m) 



(6) 



In the simplest case of a sample with (001) surface, 
where the current flows along the [100] or [010] direction, 
Eqs. (J5J reduce to 

A + B(j ■ mf, 

C{j ■m){t-m)-D(n-m). (7) 
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The magnetotransport studies presented in this work 
were performed on (001)- and (113)A-oriented samples 
with the current direction j along [110] and [332], re- 
spectively. The corresponding resistivities are: 

(001) surface and j \\ [110] 



cubic 



- A+-(B-C) + C(j-m) 2 
+ l -{C-B){n-mf, 



cubic 
r^trans 



B(j ■m)(t-m)-D(n-m). (8) 

(113) A surface and j \\ [332] 



Plo 



= A + ^(B-C) 

+ -L(126C-55)(j-m) 2 

AK 

+ —(C-B)(n-m) 2 
242 M ' 



15V2 
121 



(B-C)(j-m)(n-m), 



cubic 
r trans 



= — (9B + 2C)(j-m)(t-m) 
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The (j-m) 2 terms of pf"ng C give rise to a dependence of 
the sheet resistivity on the relative orientation between 
magnetization M and current density J, commonly re- 
ferred to as AMR. Microscopically, it is explained by 
a strong spin-orbit coupling in the semiconductor va- 
lence band. Experimentally observed differences in the 
in-plane and out-of-plane AMR, often defined as 

;0™ II 3) - Pion g (m || t) 



AMR ip = 



AMR op = 



Plongjm || J) - pip 

Plong(»™ || t) 
P\on S (m || j) - p\ong{m 



n) 



, • (io) 

have been ascribed to biaxial strain in the layer ££££1 
According to Eqs. (JSJl and ©, however, such differences 
may be expected even in the case of perfect cubic symme- 
try due to the (n • m) 2 terms of p™,^ . As will be shown 
below, a strain-induced tetragonal distortion leads to fur- 
ther (n • m) 2 terms, additionally affecting the difference 
between AMR ip and AMR op . The PHE, represented by 
the (j ■ m)(t ■ m) terms of ptrans, is closely related to 
the AMR and describes the appearance of a transverse 
voltage in the presence of an in-plane magnetic field, or 
more precisely, of an in-plane magnetization. From the 
summation terms in Eqs. (JSJ) it becomes clear that pi on g 
and ptrans, and thus AMR and PHE, not only depend 
on the relative orientation between m and j, but also on 
the orientations of m, j, and t with respect to the crystal 
axes. 

So far, quantitative studies on the angular dependences 
of the AMR and the PHE in (Ga,Mn)As were restricted 
to in-plane configurations with n ■ m = and the discus- 
sions were based on the well-known expressions^^! 

_) COS 2 , 
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Ptrans 
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(PI 



(pii — p±) sin fa cos ( 



(11) 

with <j)j denoting the angle between j and M. These ex- 
pressions, however, only hold for polycrystalline films, 
whereas (Ga,Mn)As layers are normally of high crys- 
talline quality with a uniform crystallographic orienta- 
tion of the layer. Accordingly, the expressions given in 
Eqs. are incompatible with Eqs. pif) and cannot 

be brought into agreement by simply setting n ■ m = 0. 
In fact, Eqs. ifTTl) result from Eqs. by averaging the 
summation terms over all possible crystal orientations in 
spaced with M lying in the plane spanned by j and t 

YJZ™^ = \ [2(j • ™) 2 + 1] , 



These equations are formally identical to Eqs. I|ll|) and 
apply to polycrystalline materials. A comparison be- 
tween Eqs. Ijlljl and Eqs. (|13[1 allows us to connect the 
quantities p|| and p± to the components of the galvano- 
magnetic tensors for cubic symmetry up to second order 

p\\ = A + l(W + 2C) 
5 

= Bii + r (3ami + 2an22 + 2a 232 3) , 
5 

p± = a + Ub-c) 

5 

= an + 7 (aim + 4an22 - 02323) • (14) 
5 

Thus, in general Eqs. (jl 1|> are not appropriate to describe 
the in-plane AMR and the PHE in single-crystalline 
(Ga,Mn)As layers. Only in the limiting case where 
B = C, Eqs. GJ-© simplify to Eqs. (JTTJ- 

Analyzing the angle-dependent magnetotransport data 
presented in Section IV, it turns out that additional 
terms proportional to (n ■ m) 2 have to be introduced in 
the expressions of p™ng C 1° achieve a satisfactory descrip- 
tion of the experimental results^ They are supposed to 
originate from a distortion of the crystal lattice due to 
compressive strain in the (Ga,Mn)As layer. To account 
for such strain-induced effects in a correct way, we ex- 
tend our model to a tetragonal distortion of the cubic 
lattice along the [001] direction. As will be shown below, 
our FMR and magnetotransport data suggest that this 
applies not only for the (001)- but also for the (113)A- 
oriented sample. 

In the case of a tetragonal lattice distortion along [001] , 
the symmetry reduces to D2d and the series expansion of 
p in Eq. @ yields further contributions which can be 
subsumed into an additional term Ap. The resistivity 
tensor then reads as 



with 
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-(J ■ m)(t ■ m). 
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Inserting the averaged terms into Eqs. © and using the 



relations j m — 
holds for n ■ m 
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and t ■ m 

0), we obtain 

1 ,„ 1 



sin <j)j (the latter only 



A + - [B - O + - (2B + 3(7) cos 2 fa 

O 

- (2B + 3C) sin <f>j cos fa . 



(12) The additional resistivity parameters are given by 

a = 033 — an + 01133 — an22, 
bi = CI3311 — aii22, 

&2 = 03333 — Sun — 01133 + <XH22> 
c = ai212 — 02323, 

d = a 3 i2 - ai23- (17) 

(13) Accordingly, Eqs. Q-JSJ have to be rewritten as: 
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( 001 ) surface and j 



[100] 

-\- B(j ■ m) 2 + bx(n ■ m) 2 
(C + c)(j ■ m)(t ■ m) 
- (D + d)(n-m). (18) 

(001) surface and j \\ [110] 
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(113) A surface and j \\ [332] 
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(19) 



tctra 
/^lone 



_9_ 

22 



(B-C-c) 



tctra 
r^trans 



A + 

1 

121 
9 

242 
1~21 

^(9B + 2C)(j-m)(t-m) 



11 

(126C - 5B + b + 90c) (j ■ m) 2 
(5C - bB + b + 13c) (n-m) 2 

(5B-5C -b + 9c) (j • m)(n ■ m) 



3V2 

U" 

D- 
3 a/2 



(B- 

ii 



rf(j • m), 



C) (t • m)(n • m) 
(n • m) 



(20) 



where & = 18fei + 46 2 - For perfect cubic symmetry the 
parameters a, b\, 62, c, and <i vanish and Eqs. (|18|I - H20(I 
reduce to Eqs. Q - ©- It should be emphasized that the 
expressions for pi ong and ptrans derived above generally 
apply to ferromagnets of cubic or tetragonal symmetry, 
provided that the angular dependence of the resistivity 
tensor is exclusively determined by the direction cosines 
of the magnetization. Effects correlated with the magni- 
tude B of the magnetic induction B, such as the negative 
magnetoresistance, can be easily taken into account by 
considering -B-dependent resistivity parameters. 



B. Magnetic anisotropy 

The pronounced MA in (Ga,Mn) As is associated with a 
density of the free enthalpy 28 G being highly anisotropic 
with respect to the orientation of M. The direction of 
M, i.e., the vector m which enters the equations for piong 
and ptrans given above, aligns in such a way that G takes 
its minimum. In addition to the single-domain model, 
we assume that the magnitude M of the magnetization is 
nearly constant under the given experimental conditions 
while its orientation m is strongly affected by the applied 



magnetic field H . Instead of G we therefore consider the 
normalized quantity Gm = G/M, allowing for a more 
concise description of the MA. For a biaxially strained 
(Ga,Mn)As film grown on GaAs(OOl) substrate, it can 
be written as^ 



G 



T = -LioH-m + B cll (m 4 x + m 4 y )+B c± i 
+ B 0Q1 (n ■ m) 2 + B llQ (t ■ m) 2 , 



(21) 



withn || [001] andt || [110]. The terms refer, respectively, 
to the Zeeman energy, to the cubic anisotropy under 
tetragonal distortion, to an effective uniaxial anisotropy 
perpendicular to the film including demagnetization and 
magnetoelastic effects, and to a uniaxial in-plane contri- 
bution whose origin is still under discussion ,1^22^ The 
anisotropy parameters Bi introduced in Eq. I|21|l are in 
by the aniso 
M and|H 

-^ -ff4j|/4, B cl _ = -^ #4_l/4, B- 1W = - f u H 2 \\/2, and 
-3ooi = Mo47rM e ///2. Note that by using the trivial iden- 
tity 



SI units. Expressed by the anisotropy fields Hi and 
4:irM e ff used in Refs. and l3ll they read as B c \\ 



in -' = (j ■ m) 2 + (t ■ m) 2 + (n ■ m) 2 = 1, 



(22) 



Eq. I|21|) can be easily converted to a completely equiv- 
alent expression where the in-plane contribution along 
[110] is formally replaced by a contribution along [110]. 
The only consequence is a redefinition of the anisotropy 
parameters and the addition of a constant term which 
does not alter the physical information provided by Gm- 
In Fig. [3 the free enthalpy is visualized by a 3D plot, cal- 
culated for a weak magnetic field /j,qH — 0.15 T and a set 
of anisotropy parameters with arbitrarily chosen values 
B c || = B c i_ = -0.1 T, £001 = 0.15 T, and B I10 = -0.05 T. 
The direction of M was calculated, as throughout the 
present work, by numerically minimizing Gm with re- 
spect to 9 and ip (see Fig. H). In doing so, we are able 
to trace the motion of Af, starting from a given orien- 
tation, while sweeping or rotating H . Figure [3] shows 
as an example the simulated polar (0^,0) and azimuth 
(ip H ,(p) angles of H (dashed line) and M (solid line), 
respectively, while H is rotated within the (111) plane. 
For the simulation the same field strength and the same 
anisotropy parameters have been chosen as for the 3D 
plot in Fig. [3 While H smoothly rotates within the 
(111) plane, M remains very close to the (001) plane 
(6 pa 90°) and undergoes sudden jumps in if whenever 
the minimum of Gm discontiuously changes its position. 

In the case of the (Ga,Mn)As films grown on 
GaAs(113)A substrates, the best fits to the experimental 
data (see Section IV B) are achieved for a normalized 
free-enthalpy density of the formii 



U M — 



—HoH ■ rn + B c \\ (m x + m y 
B 113 {n ■ m) 2 + B lw (t ■ mf 
Booim 2 . 



B c j_m z 



(23) 



The first five terms correspond to those already presented 
in Eq. (|21|) . The sixth term, which has to be additionally 
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[001] 



[010] 




FIG. 2: Gm as a function of m, calculated for a given mag- 
netic field /j-oH = 0.15 T and a set of anisotropy parameters 
with arbitrarily chosen values B c u = B c ± = -0.1 T, Booi = 
0.15 T, and S I10 = -0.05 T. The equilibrium position of M 
is determined by the minimum of Gm- 



150 



« 100 
T3 



a* 50 



I 1 1 1 1 I 1 1 1 
H in (111) plane 




100 200 300 

Angle of rotation a (degree) 



300 
200 8 
100 ^ 


-100 



9- 




4^ 








/ / A 



FIG. 4: The angle-dependent magnetotransport measure- 
ments were carried out for H rotated within (a) the layer 
plane, (b) a plane perpendicular to the current direction j, 
and (c) a plane spanned by j and the normal vector n. 



u>h f is fulfilled at the resonance field H res which strongly 
depends on the orientation of H due to MA. The reso- 
nance condition is given by^ 



7 



d 2 G M d 2 G 



hi 



d6 2 dip 2 



fd 2 G M \ 
\ d9dip J 



(24) 

where 7 = .9/15/1 1 denotes the gyromagnetic ratio, g 
the g-factor, [Ib the Bohr magneton, and H the Planck 
constant. At any given direction of if, the resonance 
field is obtained by solving Eq. (|24|l at the equilibrium 
position of M, i.e., for dGu/dip = and dGM/d9 = 0. 
The anisotropy parameters can then be derived from a 
fit to the measured H res recorded as a function of field 
orientation. In the present work, most of the calculations 
were carried out numerically. 



FIG. 3: Simulated polar (left axis) and azimuth (right axis) 
angles of the magnetic field H (dashed lines) and the mag- 
netization M (solid lines) for H rotated in the (111) plane. 
The same field strength and anisotropy parameters have been 
used as in Fig. |21 

introduced to obtain an optimal agreement between ex- 
periment and theory, is an inclined uniaxial contribution 
along [001], i.e., neither parallel nor perpendicular to the 
fihn>i& We attribute it to a lattice distortion along the 
[001] direction. 

C. Ferromagnetic resonance 

A highly efficient and widely used tool to study MA is 
ferromagentic resonance spectroscopy^ Most recently, a 
detailed review on FMR in (Ga,Mn)As has been given by 
Liu and Furdyna^i In the FMR experiments, the total 
magnetic moment of the Mn-ion/hole spin complex and 
thus the magnetization M precesses around its equilib- 
rium position (which in general is not identical to the 
orientation of H) at the Larmor frequency ljl- Sweep- 
ing the magnitude of the magnetic field H at a fixed 
microwave frequency u>h f , the resonance condition lul — 



IV. RESULTS AND DISCUSSION 

The longitudinal and transverse resistivities of the 
(001) and (113)A (Ga,Mn)As layers were measured as 
a function of the magnetic field orientation at fixed field 
strengths [IqH = 0.1, 0.25, and 0.7 T. In order to probe 
the anisotropy in all three directions in space, the applied 
magnetic field H was rotated within three different crys- 
tallographic planes perpendicular to n, j, and t, respec- 
tively, as shown in Fig. In the following, the measured 
angular dependences of piong and ptrans are discussed and 
values for the resistivity and anisotropy parameters are 
derived by fits to the experimental data using the the- 
oretical formalism presented in Section III. The results 
of FMR measurements, carried out on the same samples, 
are presented for reference. Note, however, that it is 
not the aim of the present study to yield a detailed or 
complete set of anisotropy parameters. In fact, the work 
is meant to provide a comprehensive theoretical tool for 
the description of the resistivies in arbitrarily oriented 
(Ga,Mn)As layers and to demonstrate the potential of 
angle-dependent magnetotransport studies for the inves- 
tigation of magnetic anisotropy. 
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(a) Angle of rotation (degree) (b) Angle of rotation (degree) (C) Angle of rotation (degree) 

FIG. 5: Angle-dependent resistivities pi on g (circles) and ptrans (dots) of the (001) (Ga,Mn) As sample at 4.2 K. The measurements 
were carried out at fixed field strengths of ^j,qH = 0.1, 0.25, and 0.7 T with H rotated in (a) the (001), (b) the (110), and (c) 
the (110) plane. The solid lines represent fits to the experimental data using Eqs. 1)19^ and one set of resistivity and anisotropy 
parameters. The dashed lines at 0.7 T simulate the limiting case where M perfectly aligns with H . In (a) the dashed lines 
completely coincide with the solid lines. 
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from a fit to the angle-dependent resistivities shown in Fig. |S] 



A. (001) orientation 

The anisotropy of the (OOl)-oriented sample was 
probed by rotating H within the (001), (110), and (110) 
planes. The corresponding angular dependences of piong 



and 



Ptr 



measured with the current direction along 



[110], are shown in Fig. [SJ At 0.7 T the Zeeman en- 
ergy dominates the free enthalpy and MA only plays a 
minor role. As a consequence, M is expected to nearly 
align with the applied magnetic field and to continu- 
ously follow the motion of H . In fact, the curves of 
Piong and ptrans at 0.7 T are smooth and largely reflect 
the anisotropy of the resistivity tensor. With decreas- 
ing magnetic field the influence of the MA increases and 
the orientation of M deviates more and more from the 
field direction. Accordingly, jumps and kinks occur in 
the curves at 0.25 and 0.1 T, arising from sudden move- 
ments of M caused by discontinuous displacements of 
the minimum of the free enthalpy. The observed angu- 
lar dependences of pi on g and ptrans can be understood 
in great detail by modeling the measured curves within 
the theoretical framework presented in Section III. For 
this purpose, the resistivity and anisotropy parameters 
from Eqs. |@J|, (|17fl . and (|21|l were determined by an iter- 
ative fit procedure. Starting with an initial guess for the 
anisotropy parameters, the resistivity parameters were 
obtained by fitting Eqs. 1|19|) to the experimental data 
recorded at 0.7 T. Then the anisotropy parameters were 
modified for an optimal agreement at 0.25 and 0.1 T, 
and the whole procedure was repeated until no further 
improvement of the fit could be achieved. The unit vec- 
tor m at any given magnetic field H was calculated by 
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numerically minimizing Gf/ in Eq. (J3TJ with respect to 
the polar and azimuth angles of M . With the excep- 
tion of A, the resistivity parameters turned out to be 
field independent within the accuracy of the fit and are 
given by B = -2.3 x 1(T 4 fi cm, C + c = -1.7 x 1(T 4 
fl cm, bi = 0.9 x 1CT 4 Q cm, and D + d = -4.4 x 1CT 4 

cm. The resistivity parameter A was found to decrease 
from 7.21 x 1CT 3 Q cm at 0.1 T to 7.08 x 10~ 3 Q cm at 
0.7 T, reflecting the negative-magnetoresistance behavior 
of pi on g.^" 7 For the anisotropy parameters we obtained the 
values B c{l = -0.015 T, B cl _ = T, B 0Q1 = 0.17 T, and 
Bi w — 0.002 T. The theoretical curves calculated with 
these parameters are in excellent agreement with the ex- 
periment and are drawn as solid lines in Fig. El 

Once the anisotropy parameters are known, the orien- 
tations of the easy axes can be determined by minimizing 
G® / with respect to m at zero magnetic field. The easy 
axes are found to lie within the (001) layer plane (0 — 
90°) at the azimuth angles ipi = 1.9° and cp% = 88.1° 
(see Fig. The slight deviation from the cubic [100] 
and [010] axes towards the [110] direction arises from the 
positive value of Bi w . 

Using the resistivity and anisotropy parameters given 
above, plots similar to Fig.|3|can be drawn for each con- 
figuration and field strength, revealing in great detail the 
motion of M. Figure El shows as an example the polar 
{9h, 8) and azimuth (tp H: ip) angles of H and M, respec- 
tively, plotted as a function of the angle of rotation a 
for H rotated in the (110) plane at 0.1 T [see Fig.EJb)]. 
While 6h passes through all values between 0° and 180° 
(dashed line, left axis), M remains very close to the (001) 
plane with 75° < 6 < 105° (solid line, left axis). For 
40° < a < 140° and 220° < a < 320°, where H is closer 
to the (001) plane than to the [001] axis, the azimuth 
angles (fn and (p of H and M, respectively, almost per- 
fectly coincide. When H approaches the hard [001] axis, 
however, M tends towards the easy [100] axis (azimuth 
angles 0° and 180°) or to the easy [010] axis (azimuth 
angles -90° and 90°). At a = _180° and a = 360°, where 
H exactly aligns with the [001] and [001] directions, re- 
spectively, the azimuth angle of M undergoes a sudden 
jump by 90°. 

Considerable information can also be obtained by com- 
paring the measured angular dependences with those ex- 
pected for the limiting case where M perfectly aligns 
with H. To this end, piong and ptrans were calculated 
with m replaced by the vector h = H/H in Eqs. (|19fl . 
The resulting curves are depicted by the dashed lines in 
Fig. For H rotated within the (001) plane [Fig. EI a)], 
the linear and quadratic (n-m) terms in Eqs. (|19|l vanish 
and we obtain the well-known cos 2 <pj and cos <pj sin <f>j de- 
pendences of piong and ptrans, respectively. At 0.7 T the 
dashed curves coincide with the solid curves, meaning 
that for H rotated within the layer plane the magnetiza- 
tion almost perfectly follows the motion of the magnetic 
field. For lower fields, M remains in the layer plane 
since [001] is a hard axis, but it increasingly deviates 
from H towards the easy [100] and [010] axes. At 0.1 T 



it abruptly switches whenever H approaches the some- 
what harder [110] and [110] axes, leading to the kinks 
observed for pi ong . 

The rotation of H within a plane perpendicular to 
the layer is accompanied by significant differences in 
the orientations of H and M, even for 0.7 T. This is 
clearly demonstrated in Figs. Etb) and (c), where the 
dashed curves represent the (n • m) 2 = cos 2 </>„ and 
(n • m) = cos0„ dependences of /9i on g and ptrans in 
Eqs. ljT§|l . respectively, with <f> n denoting the angle be- 
tween m and n. At 0.7 T, H and M coincide whenever 
H is orientated parallel or perpendicular to the layer 
plane. At lower fields this is no longer true and M avoids 
the perpendicular direction by tending towards the easy 
[100] and [010] axes (see Fig. EJ). Accordingly, the dif- 
ferences between the minimum and maximum values of 
Piong and ptrans are drastically reduced at 0.25 and 0.1 T. 

It should be emphasized that the magnctotransport 
measurements give clear evidence for the tetragonal dis- 
tortion of the (Ga,Mn)As layer: First, the parameters 
B c \\ and B c ±, representing the in-plane and out-of-plane 
contributions to the cubic-anisotropy term in G™ 1 ' s i§" 
nificantly differ. Second, Eqs. JSJ, which have been de- 
rived for the case of perfect cubic symmetry, correctly re- 
produce the measured amplitudes of piong for the in-plane 
configuration in Fig. Eta), but do not so for the two out- 
of-plane configurations in Figs.Etb) and (c). Moreover, 
it is worth noting that the fits presented in Fig. El could 
be improved even further by taking into account higher- 
order terms in the series expansion of the resistivity ten- 
sor in Eq. However, since the agreement achieved 
in second order is more than satisfactory and since the 
mathematical expressions for the resistivities would be- 
come much more complicated, higher-order terms have 
not been considered in the present study. 

The results of FMR measurements are presented in 
Fig. [3 It shows the measured and simulated angular 
dependences of the resonance field H les for H rotated 
within the (001), (110), and (110) planes. The dashed 
lines, reproducing only roughly the experimental curves, 
were numerically calculated using Eq. (|21[l. Eq. I|24fl. g — 
2.0, and the anisotropy parameters derived from the mag- 
netotransport measurements. The agreement between 
experiment and theory is significantly improved using g 
= 1.9 and the slightly higher values B c]l = -0.02 T, B c± 
= T, B 001 = 0.24 T, and B lw = 0.002 T, which were 
obtained by a least squares fit based on Eqs. I|21|) and 
(|24|l . The values agree within 30% with the anisotropy 
parameters determined from magnetotransport. The rea- 
son for the remaining difference between the two sets of 
parameters is not yet clear. Inevitable sample heating up 
to 150°C for less than 30 min during the Hall-bar prepa- 
ration as well as the different lateral sizes of our samples 
(shape anisotropy) are not expected to account for the 
observed variation. 
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B. (113)A orientation 

As already mentioned in Section III B, our experimen- 
tal data suggest the existence of a lattice distortion along 
[001] even in the (Ga,Mn)As film grown on GaAs(113)A 
substrate. This is demonstrated in Fig. |SJ which shows 
the measured and simulated angular dependences of the 
FMR field iJ rcs for H rotated within the (332), (110), 
and (113) planes. The solid and dashed lines depict the 
results of least squares fits using Eqs. (|23fl and 124fl with 
and without considering the uniaxial term i?ooi TO z7 re- 
spectively. Even though no perfect simulation of the 
measured curves could be achieved, the solid curve is 
much closer to the experimental data than the dashed 
one. Optimal agreement is obtained for g = 2.0 and 
the values B c \\ = B cl _ = -0.046 T, B 113 = 0.032 T, and 
-Sooi = 0.053 T of the anisotropy parameters (solid line 
in Fig.[SJ|. Similar to the case of the (001) layer, the uni- 
axial in-plane contribution is almost negligible with Bj 10 
= -0.005 T. 

The magnetotransport data, measured with the cur- 
rent direction along [332], are depicted in Fig. The 
figure shows the experimental and simulated angular de- 
pendences of piong and ptrans for H rotated within the 
(113), (332), and (110) planes. Using the fit procedure 
described in Section IV A, the values of the resistivity 
parameters in Eqs. H20f) are obtained as B — —8.1 x 10~ 4 
ft cm, C = -8.7 x 10~ 4 fl cm, b = 7.2 x 10~ 3 fl cm, 
c = 5.1 x 10~ 4 n cm, D = -2.1 x 10~ 3 ft cm, and d 
— 2.1 x 10~ 4 ft cm. The sum A + 2a/ll decreases from 
21.61 x 10~ 3 ft cm at 0.1 T to 20.79 x 10" 3 ft cm at 0.7 T. 
For the anisotropy parameters we obtained the values 
B c]] = -0.046 T, B c± = -0.03 T, B 113 = 0.018 T, B 001 = 
0.02 T, and B ll0 = -0.008 T. Similar to the (001) sam- 
ple, they are smaller than the values determined by the 
FMR study, however, the discrepancy is less than a fac- 
tor of about two. The calculated curves, represented by 
the solid lines in Fig. are in excellent agreement with 
the measured data. Again, the dashed curves simulate 
the case of a magnetization which perfectly aligns with 
H. A comparison between the dashed and solid curves 
reveals that the magnetization at 0.7 T almost perfectly 
follows the motion of the magnetic field in the two out-of- 
plane configurations [Figs.|5{b) and (c)], whereas for H 
rotated within the layer plane [Fig.^a)] M significantly 
deviates from H. According to the model calculations, 
the latter behavior arises from the cubic terms and the 
uniaxial [001] contribution in G 1 ^ 3 , resulting in a deflec- 
tion of M towards the (001) plane. The asymmetry of 
Ptrans in Fig. Ela) partly results from this deflection and 
partly from the last term in Eqs. (|20fl which originates 
from the tetragonal distortion. A similar asymmetry has 
been observed by Muduli et al^» in FeaSi films grown 
on GaAs(113)A substrates. There, the asymmetry has 
been explained by third-order terms in ptrans- We can- 
not rule out that in the (113)A (Ga,Mn)As sample un- 
der study higher-order terms contribute to the asymme- 
try, too. However, since the present model, including 
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terms up to second order, fully accounts for the observed 
angular dependence of pi on g and Ptransj it has not been 
considered as mandatory to include them. 

Using the anisotropy parameters obtained from the 
curve fits, the orientations of the easy axes were deter- 
mined by minimizing G\] 3 with respect to m at zero 
magnetic field. We find the easy axes in the (113) A layer 
at the angles (see Fig. Q} 6 X = 92°, ip x = -3.4°, and 9 2 
= 92°, tp 2 = 93.4°, i.e., very close to the [100] and [010] 
axes, in qualitative agreement with the results presented 
in Refs.Elandlil 



SUMMARY 



symmetry The expressions, applicable to (Ga,Mn)As 
layers with arbitrary surface index, were used to quanti- 
tatively model the angular dependences of the resistiv- 
ities, measured in (001) and (113)A (Ga,Mn)As films 
as a function of magnetic field orientation. Whereas 
the curves at 0.7 T largely reflect the anisotropy of 
the resistivity tensor, the curves at 0.25 and 0.1 T are 
strongly affected by magnetic anisotropy, allowing access 
to anisotropy parameters. The magnetotransport data 
and comparative ferromagnetic resonance studies reveal 
an inclined uniaxial anisotropy along [001] in the (113)A- 
oriented (Ga,Mn)As layers in addition to the usual in- 
and out-of-plane contributions known from (001) layers. 



A series expansion of the resistivity tensor with respect 
to the magnetization components yields general expres- 
sions for the longitudinal and transverse resistivites in 
single-crystalline ferromagnets with cubic and tetragonal 
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